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Abstract 

We present Lagrangian one-particle statistics from the Ris0 PTV experiment of a turbulent flow. 
We estimate the Lagrangian Kolmogorov constant Co and find that it is affected by the large scale 
inhomogeneities of the flow. The pdf of temporal velocity increments are highly non-Gaussian for 
small times which we interpret as a consequence of intermittency. Using Extended Self-Similarity we 
manage to quantify the intermittency and find that the deviations from Kolmogorov 1941 similarity 
scaling is larger in the Lagrangian framework than in the Eulerian. Through the multifractal model 
we calculate the multifractal dimension spectrum. 



I. INTRODUCTION 



In the present contribution we present experimental results on Lagrangian one-particle 
statistics from an experiment with the Ris0 Particle Tracking Velocimetry (PTV) setup. We 
focus on small-scale statistics in a turbulent flow: the statistic is analyzed with Extended 

n 

Self- Similarity (ESS) [2| and the results are presented in the spirit of the multifractal model 
of turbulence The use of ESS is discussed together with the multifractal model in a 
finite Reynolds number flow like the present. 

We have performed a Particle Tracking Velocimetry (PTV) experiment in an intermediate 
Reynolds number turbulent flow. The flow has earlier been reported in |3,U,Q] although we 
use data from a recording with a slightly lower Reynolds number. PTV is an experimental 



20 



24, 



method suitable for obtaining Lagrangian statistics in turbulent flows |9|, 
Lagrangian trajectories of fluid particles in water are obtained by tracking neutrally buoyant 
particles in space and time. The flow is generated by eight rotating propellers, which change 
their rotational direction in fixed intervals in order to suppress a mean flow, placed in the 
corners of a tank with dimensions 32 x 32 x 50cm 3 (see Fig^). The data acquisition system 





FIG. 1: Experimental setup 
consists of four commercial CCD cameras with a maximum frame rate of 50Hz at 1000 x 1000 
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TABLE I: Turbulence characteristics: e is the mean kinetic energy dissipation, 77 = (^ 3 /e) 1 / 4 is the Kol- 
mogorov length scale with the kinematic viscosity v = 1 of water. t v = {v/e) 1 / 2 is the Kolmogorov time 
scale and cr^ = | (cr 2 .^ + cr\ + a\ z ) is the standard deviation of velocity. The integral length scale is de- 
fined as L = a 3 /e while Tg is the eddy turnover time Te = L/a u . The Reynolds number is defined as 
Re x = Vl5(L/ V ) 2 ^ 3 . 



pixels. The measuring volume covers roughly (12cm) 3 . We use polystyrene particles with 
size ~ 400/zm and density very close to one. We follow 0(1000) particles at each time step 
with a position accuracy of 0.05 pixels corresponding to less than 10/im. 

The Stokes number, 77 / t v (tj denotes the inertial relaxation time for the particle to the 
flow while t v is the Kolmogorov time) is much less than one and the particles can therefore 
be treated as passive tracers in the flow. The particles are illuminated by a 250W flash 
lamp. 

The mathematical algorithms for translating two dimensional image coordinates from the 
four camera chips into a full set of three dimensional trajectories in time involve several cru- 
cial steps: fitting gaussian profiles to the 2d images, stereo matching (line of sight crossings) 



with a two media (water-air) optical model and construction of 3d tra 
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ectories in time by 



using the kinematic principle of minimum change in acceleration 

The flow characteristics are presented in Table |U The mean flow is axisymmetric with 
a significant vertical straining on the largest scales and we did not find any significant 
differences from the flow reported in ^, ^| , where properties of the mean flow can be found. 

Here we look at a sub-volume of the full measuring volume. Only particles which we 
can observe within a ball of radius 50mm is considered and the turbulence characteristics 
given in Table U] are thus only determined from particles inside this ball. The ball is centered 
approximately in the center of the tank where the velocity standard deviation a u has a global 
minimum. Inside the ball the particles are uniformly distributed. With r v = 0.09s and a 
recording frequency at 50Hz the temporal resolution is ~ 4frames/r J? . 

The database is the largest we have compiled and it consists of ~ 10 6 individual tra- 
jectories with an average length of ~ 8^, a standard deviation of ~ 13r^ and the longest 
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tracks we find are ~ 150t v . The number of tracks was an important requirement since the 
calculation of high order moments is considered important. 

Throughout the paper we will denote the Lagrangian velocity along a particle trajectory 
for v(i) and the Eulerian velocity in a fixed frame of reference for u(x, t). 



II. FINITE VOLUME MEASUREMENTS 



A nice property of the Eulerian velocity statistic is that it is stationary in time in the 
present experiment. This is not the case for the Lagrangian statistics. The non-stationarity 
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FIG. 2: (v 2 (t + to)). The average is taken over all particles which were observed inside B at both time to 
and to + 1. 



is showed in Figure 121 where (v 2 ) is observed to decrease over time. This reflects the finite 
measuring volume and the non- uniform forcing in space in our experiment: the particles only 
gain kinetic energy close to the propellers. During their subsequent motion the particles loose 
kinetic energy until they again come close to the propellers which are constantly spinning. 
Looking at a finite measuring volume away from the propellers, there will therefore be a flux 
of kinetic energy into the volume. Inside the volume the kinetic energy is dissipated and 
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hence we have at the entry of the volume 

\^ = w 

where we recognize the mean kinetic energy dissipation e. From Figure El we find e = 
124mm 2 /s 2 . This number is close to the number obtained from the second order Eulerian 
Structure Function e = 132mm 2 /s 2 . We take the difference as the uncertainty in estimating 
e. Eqn. [T]can also be derived directly from the Navier-Stokes equation by assuming global 
homogeneity. 

The Lagrangian second-order structure function is defined as 

S 2 (r) = ([v(t + r)-v(t)] 2 ), (2) 

where v(t) is here the velocity component along a fluid trajectory. Similar the Lagrangian 
co-variance function is defined as 

R L (r) = (v(t)v(t + r)). (3) 

The non-stationarity of (v 2 ) means that 

Sl(t) = (v 2 (t)} + (it 2 ) - 2R L (t) < 2((u 2 } - R L (t)), (4) 

where we have used that the Lagrangian velocity on the boundary of the measuring volume 
B equals the Eulerian velocity and therefore (v 2 (t = 0)) = (u 2 ). S\{t) is plotted in Figure El 
for all three velocity components. It is clear that for long times S\ does not approach 2{u 2 ) 
in agreement with eqn. |3J 

A common interpretation of the finite volume influence on Lagrangian statistics is that 
the particles we can observe for long times are relatively slow ones or particles which are 
trapped in high intensity vortices (see later). Here we emphasize the equivalence with the 
energy argument of decaying turbulence described above. Particles which can be observed 
for long times are slow because it is long time ago they gained kinetic energy at the forcing 
site. 

In Direct numerical simulations (DNS) forcing occurs in wave-number space on the lowest 
wave-numbers. We therefore have d(v 2 )/dt = and consequently Lagrangian stationarity. 
Most physical flows encountered in nature will, however, be Lagrangian non-stationary. 
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FIG. 3: Second order Lagrangian structure function for the three coordinates of v(t). x : green (radial 
component), y : red (radial component) and z : blue (vertical component). The horizontal lines is the 
Eulerian velocity variance (u 2 ). 



III. ANISOTROPY AND INERTIAL RANGE SCALING 

The linear dependence of Re\ on Tl/t^ implies that a very high Reynolds number is 
needed in order to obtain a clear Lagrangian inertial range. Yeung concluded, based 
on extrapolations from Eulerian fields in DNS, that at least Re\ ~ 600 — 700 was needed. 
Experimental flows at Re\ = 1000 ^| and Re\ = 815 do, however, not show a very 
pronounced inertial range defined as a linear regime in the second-order structure function 
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In the inertial range t v < t < Tl, K41 similarity theory predicts 

S 2 L>ij (T) = CoerSij, (5) 

where the Lagrangian Kolmogorov constant Co is supposed to be universal for infinite 
Reynolds numbers j^. Cq plays a crucial role in stochastic models j3| and has lately 
been shown to reflect anisotropy in the large-scale forcing [2^. In Figure E] we present re- 
sults of C for the three coordinates of v(t). According to eqn. C should be determined 
from a plateau in the inertial range. The parabolic form therefore reflects the almost vanish- 
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FIG. 4: Cq for the radial components (green and red) and the axisymmctric component (blue). 



ing inertial range in our experiments. The difference between radial and the axisymmetric 
component stems from the large scale anisotropy. Since Co is maximum for times around 
2 — At and therefore mainly associated with small scales the difference is a clear signature 
of small-scale anisotropy. The values of Co are 5.34 ± 0.16, 5.08 ± 0.15 and 4.09 ± 0.12 for 
the three components x, y and z respectively. 

It is interesting to see that the slight difference in the radial forcing is surviving all the way 
down. The propellers forcing the flow are attached to four rods placed in the corners of the 
tank. The reason for the radial components being different is probably small differences in 
the manual vertical placement of the propellers on the rods. The lack of small-scale isotropy 
in the current experiment should not be taken as a failure of Kolmogorov's hypothesis of 
local isotropy. For that the Reynolds number is not high enough. Other experiments at 
much higher Reynolds number do, however, all indicate that the large scale inhomogeneities 
are also present at smaller scales although with smaller amplitude 22, 28, 29^ . 

Alternatively one can calculate the lagrangian velocity spectrum $(o;) and calculate Cq 
from this. is defined as the fourier transform of the velocity co- variance function 

Rl(t) 0: 

1 f°° 

$(w) = — / dTexp(-iuT)R L (T). (6) 
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In the inertial range K41 predicts 

§ l3 {u) = f3eu-\ (7) 

with Co = 7i/3. In Figure El (a) we have plotted Ri(t) in the three directions. The radial 
components fall off exponential with e-folding times T* xp ~ 10. 7t v and Tf xp ~ 9.4-r^ while 
the vertical axisymmetric component T* xp ~ 14^. Since -Rl(t) is composed of eddies of 
all size in the flow, the energy containing scales and hence the large scale inhomogeneities 
strongly effects its form. The integration of Ri{t)/a 2 gives the Lagrangian integral time 
scale T L . We find values of T L ~ T exp . The velocity spectrum $(u;) is shown in (b). For 




FIG. 5: (a) Rl{t). (b) The straight line is the K41 prediction ~ id 2 . Color codes as in Figure^ 

The curves have been shifted vertically for clarity. 

small frequencies u the spectrum for all three components are white. This corresponds to 
uncorrelated velocities for long time lags on a track. For higher frequencies all three spectra 
turn red with slope of ~ —2 in agreement with the Kolmogorov prediction. Due to a relative 
low sampling rate (dt = 0.021s) the Nyquist frequency prevent us from studying frequencies 
higher than uj = 23.8s -1 . 

Lien and D'Asaro studied the scaling properties in a simple Lorentzian model spec- 
trum and found that with a finite Reynolds number it is easier to obtain inertial range scaling 
from the spectrum than from the structure function 5£(r). We have plotted the spectrum 
compensated with uj 2 in Figure |H1 in order to have a better look at the existence of an iner- 
tial range. For all three components a narrow inertial range is observed as a plateau. The 
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horizontal lines are used for estimating Co- We find values equal to 4.91 ± 0.15, 4.79 ± 0.14 
and 4.07 ± 0.12 for the three components respectively. These values are smaller and a bit 
more isotropic than those calculated from the structure functions. This is in contrast to the 

n 

arguments by Lien and D'Asaro [12J C$ who state that they should be larger. 
10000 




FIG. 6: Compensated velocity spectrum uj 2 ^(uj). Color code as in Figured The curves have been shifted 
vertically for clarity. One can therefore not determine the magnitude of (jJ 2 &(uj) from the different curves. 
The horizontal lines are the levels from which Co is calculated. 



IV. SMALL-SCALE INTERMITTENCY 

From the study of the lower moments we proceed to higher order moments describing the 
most extreme events. 

The pdfs of temporal velocity increments Sv(t) = v(t + r) — v(t) are shown in Figure [7| 
for different time lags r. All three components are shown. The three components show the 
same over all behavior: for large time lags the distributions are Gaussian while they for 
small time lags have fat tails. The curves corresponding to the smallest time lags have a flat 
plateau at 5v ~ 0. This is a binning artifact and does therefore not represent any physical 
trend in the data. The non-Gaussianity for small times becomes more clear by looking at 
the flatness. The flatness of the distributions is defined as 
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FIG. 7: pdf of velocity increments Sv(t) for times (decreasing downwards) t = 
1.0t v , 1.7t, ; , 3.6t v , 6Trj, 12t v , 2At v and 48t v . The curves have been shifted vertically for clarity. Color 
coding as in Figure 01 The black curve is a Gaussian. 



F(r) 



(8) 



{5v 2 (t)) 2 

and is shown in Figure |H1 F is monotonically decreasing for all three components and reaches 
a Gaussian level at time lags: r ~ 40^, which is substantial larger than T^. We do not at 
present have any explanation for this. 

The results presented in Figure |7] and |H] are strong evidence of Lagrangian intermittency, 
i.e. non-Gaussian behavior of the smallest temporal scales in the flow. These results agree 

ri 

with observations by Mordant et al. |l9| and Direct Numerical Simulations (DNS) by Biferale 
et al. p]. 

Our findings suggest that intermittency can be studied in flows with a moderate Reynolds 
number of order (9(100). The only necessary condition seems to be the size of the ensemble: 
a large number of particles is needed to observe rare events. 
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FIG. 8: The flatness of Sv(t). The color coding as in Figure^] The horizontal line F = 3 is the Gaussian 
prediction. 

A. Higher order structure functions and ESS 

Before we look at the higher order moments we check for convergence of these. In Figure El 
we show 8v u (t)p(5v(t)) for n = 4,6,8, 10. The time lag in all four plot is r = 2.1t v . For 
n < 8 we observe convergence. For n = 8 we start to get into trouble, but it seems like we 
have captured most of the signal - at least for the radial components (red and green curves) . 
In an incompressible flow (5u(t)) = 0: the non-zero skewness observed in all the curves is 
therefore an artifact of sampling errors and / or tracking of particles. This is an issue which 
has to be resolved before more dramatic conclusions can be made. 

K41 similarity theory predicts in the inertial range for the p order structure function: 



Intermittency can be defined as the departure from K41 similarity scaling. This means that 
eqn. |U]can be replaced by a more general form taking intermittency into account: 



S p L (r) = (5v p {t)) ~e p/2 t p ' 2 



(9) 




(10) 



where is commonly named the Lagrangian anomalous scaling exponent. Only recently 
it has been possible to measure and hence quantitatively describe the extreme dynamics 
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FIG. 9: 6v n p(6v) for n = 4, 6, 8, 10. The time in all four plot is r = 2.1t v . 
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. The data 



present in the fat tails of the distribution of 8v{r) for r — > 
presented here is therefore merely a verification of already obtained results. 

In Figure El (a) structure functions S^(t) of order p = 2,4, 5, 6, 8 are shown as a function 
of t. Power laws have been fitted to each function in the region 2r v <t< 4t v corresponding 
to the maxima of S p L {t)/t/e. The fits are not convincing. First of all, the inertial range 
is too narrow and we therefore can not expect any universal scaling. Secondly, and less 
importantly, we know that the small scales are affected by the large-scale inhomogeneities. 

A popular way of looking at scaling exponent is instead to measure ratios of scaling 
exponents. This method is called Extended self-similarity (ESS) and was introduced by 
Benzi et al. |2j . The method was introduced in the Eulerian frame but can be transferred to 
the Lagrangian frame if we assume that = 1 following K41 similarity theory. The crucial 
step is to treat all velocity increments as positive. This affects the odd-numbered structure 
functions. We therefore define 



,-L,ESS 



(\Sv(t)\ p ) ~ (6v 2 (t))^' "". (11) 

In Figure ITU1 (b) S P L ESS (t) is shown as a function of S^it). The scaling is now much better, 
which explains the wide popularity of the method. The different scaling exponents are 
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FIG. 10: (a) S£(t) as a function of t/r v for p = 2,4,6,8 increasing upwards, (b) Extended self-similarity: 
S£ ES g(t) as a function of S\{t) for p = —1, 1, 2, 9. In both panels data from the radial ^-component are 
used. 



p 


-1 


1 


2 


3 


4 








0.98 ±0.07 




1.40 ±0.08 


A L,ESS 


-0.62 ±0.07 


0.59 ± 0.02 


1 


1.27 ±0.03 


1.43 ± 0.06 


P 


5 


6 


7 


8 


9 






1.56 ±0.10 




1.66 ±0.17 




j-L,ESS 
^P 


1.53 ±0.05 
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1.75 ±0.19 



TABLE II: Lagrangian anomalous scaling exponents 



printed in Table |TI] and plotted in Figure ^2 for the radial components. The error bars 
represent small deviations between the two radial components as well as an error estimated 
from fitting the straight lines in Figure E3 The errors increase with p and are significantly 
larger for (r compared to the ESS approach (p ,ESS . 

Some remarks about ESS should be made at this point. In the original paper Benzi et al. 
Q] argued, based on experimental evidence of \(5 r u 3 (r))\ ~ (\5 r u(r)\ 3 ), that (\5 r u(r)\ p ) = 
B p {\S r u(r) l 3 )^ . As also emphasized in the paper this is not a rigorous result which can be 
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FIG. 11: Lagrangian anomalous scaling exponents Cp (blue curve) and the Extended self-similarity anoma- 
lous scaling exponents (p' ESS (red curve). The straight line is the K41 prediction. 



deduced from the Navier-Stokes equation. By plotting absolute (defined by positive velocity 
increments) structure functions vs. the third order structure function (Eulerian frame) or 
the second order structure function (Lagrangian frame), an extended scaling range can be 
observed because undulations in the structure functions are correlated and hence disappear 
when plotted against each other. ESS is widely used and gives seemingly universal scaling 
exponents for flows in a wide range of Reynolds numbers. As pointed out by Arneodo et al. [l| 
no consensus besides the observed facts exists about the interpretation or even significance 
of ESS. Whether the observed scaling in ESS is the signature of hidden physical laws is 
speculated. In the Lagrangian frame an additional problem arise. As already mentioned 
K41 predicts linear dependence of time scale for the Lagrangian second order structure 
function and hence — !• This is motivated by the scaling in the Eulerian frame and 
specifically from the four-fifth law. A similar exact result does not exist for the Lagrangian 
structure functions. So all in all, one could state that it is a wonder that it works! 

The va. ues in TabMJare in excellent „t with results obtained b y Xu et al. Q 
and Mordant et al. |l sj] . The values by Biferale et al. are somehow higher and was by 
Xu et al. (33] explained as a different choice of inertial range. 
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B. The multifractal framework 



The multifractal model of turbulence was introduced by Parisi and Frisch 231 in the 

n 

Eulerian frame after an early attempt by Mandelbrot [15J who used multifractal measures 
to characterize the spatial structure of dissipation in turbulence. 

The multifractal model is phenomenological and has been able to successfully predict the 
corrections to K41 similarity scaling for high order moments of spatial velocity increments 
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Borgas 8j discusses multifractals in the Lagrangian frame and introduces a bridge to 
the Eulerian framework. The literature is, however, not very rich on work on Lagrangian 
multifractals, which could have to do with the difficulties in obtaining reliable Lagrangian 
data set more than a animosity against the multifractal model. Work by Biferale et al. 
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34j have, however, shed light 



Chevillard et al. [lfj, Mordant et al. |12lllS|, Xu et al. 
on the issue of multifractals in the Lagrangian frame. 

In the Lagrangian multifractal model the flow is assumed to possess a range of scaling 
exponents h min , h max with a certain probability so that 

5v{t) ~ r h . (12) 

For each scaling exponent h there is a fractal set with a /i-dependent dimension D L {h). The 
embedding dimension is one (r G TV) and hence D L (h) < 1 for all h. The probability P^{t) 
of having an exponent h at time r is therefore proportional to 1 — D L (h). From a steepest 
descent argument one can calculate a relation between the anomalous scaling exponents 
and the fractal dimension D L (h) given by [ill ]: 

g = m£\ph+l-D L (h)). (13) 
If D L {h) is concave a Legendre transformation gives 

D L (h) = inf[p/i + l-C L ]. (14) 

p 

In Figure El we have plotted D L (h) obtained through eqn. EI First we calculated for 
both integer and non-integer values of p between p = — 1 and p = 9. The result is the red 
curve in the Figure fT2l 

The black dots are the result by Xu et al. [3^ who in a PTV experiment of Reynolds 
number Re\ = 200,690 and 815 measured D L {h) both trough P^(r) which they manage 
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FIG. 12: The Lagrangian multifractal dimension spectrum D L (h). The black dots are the result by Xu 
et al. tail. 



to measure directly and through eqn. El as we have done here. They arrived at the same 
D L (h) from both calculations putting confidence in the multifractal model for Lagrangian 
velocity increments. The agreement between their data and ours is very good. Only for 
h > 0.6 in the linear portion of D L (h) do we observe a discrepancy. 



This linear portion of D L (h) was by Xu et al. 
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explained in the following way: because 



the domain of h is finite eqn. 1131 will become a linear function after some p'. This linear 
behavior is also observed in Figure ^2 for large ps. The transition point h' happens where p' 
minimizes the right hand side of eqn. HH1 For p > p' we therefore have that = h min p + 1. 
Since only moments of the structure functions of order larger than —1 converge we have 
p' = -1 and the linear part of the curve is D L (h) = -h + 1 - C-i- Xu et al. H 

successfully 

corrected the models by Biferale et al. jsj (from a theoretical prediction by She and Leveque 
27 1) and Chevillard et al. [10] and found a remarkable match. The discrepancy in Figure 
IT2"1 therefore stems from different estimates of and the uncertainty in measuring it. 

Chevillard et al. |10J] came up with a formula for the connection between D L (h) and its 
Eulerian counterpart D E (h). The formula is 

h 



D (K) = -h + (1 + h) [D 



1 + h 



- 2 



(15) 
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From our database we have calculated the Eulerian anomalous scaling exponents from ESS 
structure functions 

S P E ,Ess(r) = (\tMr)n ~ (\5 r u(r)\ 3 )^ ESS ■ (16) 

Results are shown in Figure IT31 In (a) S E ESS (r) are plotted from p — 1, ...9. For all orders 
ESS seems to work fine. In (b) the anomalous scaling exponents C E,ESS are shown. 

The values are in perfect agreement with the theoretical model by She and Leveque (27 1. 
More interesting is the departure from the K41 prediction which is smaller than in the 
Lagrangian frame. This is interpret as Lagrangian statistics being more intermittent. 

Just like in the Lagrangian frame there is a Legendre transformation between Cf and 
D E (h): 

D E (h)=mf\ph + 3-( E ] (17) 
p y 

The only difference from eqn. El is the embedding dimension which in the Eulerian frame 
is three (r G 1Z 3 ). 
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FIG. 13: (a) ESS Eulerian structure functions ESS (r) of order p as a function of S^(r). p = 1,..,9 
increasing upwards, (b) Anomalous scaling exponent determined from (a) (red dots). The straight line is 
the K41 prediction and the green curve is the theoretical model by She and Leveque 27 1. 



From eqn. Eland E3 w e can find D L (h) from the Eulerian anomalous scaling exponent 
presented in Figure HH1 (b). The comparison is plotted in Figure H3I Again we observe 
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a discrepancy in the linear part of D%. Whether it comes from the determination of the 
anomalous scaling exponents from ESS or that there is a flaw in eqn.^Jwe can not say at the 
moment. A direct measurement of the probability in both the Eulerian and Lagrangian 
frame might give more insight into the connection between the two frames. 
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FIG. 14: The Lagrangian multifractal dimension spectrum D(h). The red curve is calculated from La- 
grangian ESS structure functions while the blue is obtained though Eulerian structure functions and 1151 

The physical interpretation of the multifractal model is not that easy. In K41 similarity 
scaling only one scaling exponent is possible, namely h = 1/3 and thus = p/3. This is 
motivated from the fact the Navier-Stokes (N-S) equation is only invariant under one scaling 
group. This group is characterized by an exponent h obtained by scaling the N-S with the 
following transformation (time, position, velocity): t, r, u i— > \ 1 ~ h t, Ar, A^u for A G 1Z + . The 
solution is h = — 1. In the limit of infinite Reynolds number the viscosity term in the N-S 
equation becomes negligible and we find that the N-S equation, is now invariant to infinitely 
many exponents h. This is one of the motivations for the multifractal model. It is, however, 
not a justification. Another important aspect of the model is the fact that when an eddy 
breaks up into smaller eddies in the Richardson picture the smaller eddies do not cover the 
same amount of space. Instead they cover only a fraction equal to 3 — D E {h). We thus 
have regions in the flow with large activity and regions with almost calm waters. In the 
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Lagrangian frame this would mean that the individual fluid particles are not free to move 
around in all directions. For example as reported by [24] and p| are particles often trapped 
by intense vortices. The large accelerations and velocity increments of these events are 
therefore of dimension lower than 3 in the Eulerian frame and lower than 1 in the temporal. 
This spiral motion of fluid particle around a fluid filament is also the fluid mechanical picture 
of intermittent events in the model by She and Leveque j^: by entraining surrounding fluid 
kinetic energy fluctuations are effectively dissipated along the axis of the filament. 

As emphasized by Borgas [^] the multifractal model does, however, not imply that the 
trajectories of fluid particles are fractal trajectories themselves. 



V. CONCLUSIONS 



We have measured Lagrangian one-particle statistics and looked at small-scale behavior. 
The finiteness of the measuring volume can be used to calculate the mean kinetic energy 
dissipation e in the flow without any further assumptions. The small scales do seem to 
be affected by the large-scale inhomogeneities present in our flow. We do not observe a 
significant inertial range but by Extended Self-Similarity we are able to extract a quantitative 
measure of the structure functions of high order. From these we calculate the Lagrangian 
anomalous scaling exponents and find excellent agreement with already published results. 

Via the multifractal model we have calculated the Lagrangian multifractal dimension 
spectrum. The spectrum is similar to the one published by Xu et al. j^] even though our 
Reynolds number is significantly lower and our mean flow is different. 

Most importantly we have shown that a high Reynolds number is not necessary to obtain 
results in the Lagrangian frame. All experiments and DNS do show the same qualitative 
features and no clear Lagrangian inertial range has been observed. Whether it is because 
current experiments are performed with too low Reynolds number or it simply do not exist 
future experiments will tell. 
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